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Abstract: 

We study a three-parameters family of solutions of the Brans-Dicke field equations. 
They are static and spherically symmetric. We find the range of parameters for which this 
solution represents a black hole different from the Schwarzschild one. We find a subfamily 
of solutions which agrees with experiments and observations in the solar system. We 
discuss some astrophysical applications and the consequences on the "no hair" theorems 
for black holes. 
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1. Introduction 



Lately there have been some renewed interest in the Brans-Dicke theory of gravita- 
tion. On one hand, it has been applied to cosmological models of the universe during the 
inflationary era to make more natural bubble percolation W. Also, it was found that in 
the low-energy regime, the theory of fundamental strings can be reduced to an effective 
Brans-Dicke one^. The subject of gravitational collapse, however, has not yet been thor- 
oughly studied. One of the outstanding results on this field is the Hawking theorem! 3 ], 
that states that the Schwarzschild metric is the only spherically symmetric solution of 
vacuum Brans- Dicke field equations. The proof of this theorem goes through the fact that 
the Brans-Dicke scalar field 4> must be constant outside the black hole and the use of the 
weak energy condition. In this paper we study a three-parameters family of solutions of 
Brans-Dicke equations which is static and spherically symmetric. We study under which 
range of the parameters we can have non-singular (at the horizon) black hole solutions. 
We are able to obtain explicit examples where the metric represents a black hole solution 
different from the Schwarzschild one: 

ds 2 = -A{ry- n dt 2 + AirY^dr 2 + r 2 A{r) n <m 2 ; A(r) = 1 - 2— ; n < -1 . (1) 



r 



where ro is an arbitrary constant and n represents a scalar hair. 

Classical scalar hairs in General Relativity Black Hole solutions have already been 
found for several coupling. These include the case of an axion with an RR coupling W, 
and a dilaton ^ and an axionlike scalar field ' 6 ' coupled to Einstein Maxwell theory. A 
conformally coupled scalar field can have a static ^, but unstable ^ solution. 

The Brans-Dicke theory ^ incorporates the Mach principle, which states that the 
phenomenon of inertia must arise from accelerations with respect to the general mass 
distribution of the universe. This theory is self-consistent, complete and for \u\ > 500 
in accord with solar system observations and experiments ^ 10 \ It is, in some sense, the 
simplest extension of General Relativity. It introduces an additional long-range scalar field 
(f> besides the metric tensor of the spacetime g^ v from which are constructed the covariant 
derivative and the curvature tensors, in the usual manner. u> is the Dicke dimensionless 
coupling constant. 

The theory is metric, i. e. the weak equivalence principle is satisfied. The matter 
couples minimally to the metric and not directly to (p. The scalar field does not exert any 
direct influence on matter, its only role is that of participating in the field equations that 
determine the geometry of the spacetime. 

The action for the Brans-Dicke theory is: 

s = J dx 4 V=g[<i>R-u(<p, a <f>' a )/<P + ^L matter } . (2) 

The variational principle gives the field equations: 
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The matter stress-energy tensor and together generate the metric. The field equation 
for is: 

*»= D * = iifb T - (4) 

In the next section we study a solution^ 9 ! of the vacuum field Eqs. (3)-(4). This is a three- 
parameters static spherically symmetric metric. We study the asymptotic behavior, the 
occurrence of singularities and event horizons. In the third section we study the special 
cases for which this metric can be of astrophysical relevance. We compute the geodesies 
equations, post Newtonian parameters, energy, period and redshift of the last stable cir- 
cular orbit, dispersion cross sections, Kruskal transformations and Hawking temperature. 
We end the paper with the discussion of the obtained results, in particular, the relevance 
of the non-Schwarzschild-like black holes found in the Brans-Dicke theory. 



2. Static Spherically Symmetric Vacuum Solutions 

The Brans-Dicke vacuum field equations can be written as: 



R a f3 = ^0,a0,/3 + , (5) 



□ = . (6) 

It is easy to show that a power generalization of the Schwarzschild metric is a solution of 
this equations I 9 ' 11 !; 

ds 2 = -A{r) m+1 dt 2 + A{r) n ~ 1 dr 2 + r 2 A{r) n dn 2 ; (7) 

dtt 2 = d-d 2 + sin 2 dip 2 ; 

Air) = 1-2— , 
r 

and the scalar field: 

0(0 = 00^)-^ ; (8) 
where m, n, 0o and ro are arbitrary constants. The coupling constant is found from: 

(m 2 + n 2 + nm + m — n) 

10 = —2 -7; . (9) 

[m + n) z 

Either from (7) or (8) we can compute the components of the Ricci tensor: 
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R 00 = (m + l)(m + n)^A(rr- n , (10) 



2 

R u = (- m 2 + nm + 3 n - m )^A(r)- 2 + 2(m + n)^A(r)- 1 , (11) 
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R 22 = - n (m + n) f^j A(r)" 1 - (m + n)^ , (12) 
R 33 = sm 2 i)R 22 , (13) 

And the curvature R is given by: 

2 

# = -2(m 2 + n 2 + mn + m-n)^|A(r)- n - 1 . (14) 

We observe that it vanishes like r -4 as r — > oo.We will study now the geometrical properties 
of the metric (7) for given values of the parameters m and n. 

To see that the metric (7) is asymptotically flat it is enough to show that the metric 
components behave in an appropriate way at large r-coordinate values, e.g., g^ v = + 
0(l/r) as r — > oo.By inspection of the coefficients, we verify that this is so. No matter 
which power of A(r), can be written as a binomial series: 

A(r)< = (l - 2 r f) q = 1 - ^ + q ( q - 1) (^) * + .... (15) 

Thus, asymptotically flatness is verified for every value of m and n. 

To study the occurrence of true singularities of the metric (7), (not coordinate sys- 
tem pathologies), it is enough for us to examine scalars formed out of the curvature. In 
particular, the scalar invariant: 

T D JDCe/3-fS 

1 — -Ka/3 7 <5-K — 

= 4r 2 (r - ro )-2(n+D r -4+2n . j ' h ( m , n ) + 4 (^) I 2 (m, n) + 67 3 (m, n) J , (16) 
where: 

7i(m, n) = 48 + 56m + 41m 2 + 10m 3 + m 4 - 56n - 34mn+ 
-20m 2 n - 2m 3 n + 29n 2 + 6mn 2 + 3m 2 n 2 - 8n 3 + n 4 , (17) 
I 2 (m,n) = -12 - 13m - 8m 2 - m 3 + 13n + Aran + 2m 2 n - 6n 2 + n 3 , (18) 
h{m,n) = (m+l) 2 + (n-l) 2 . (19) 

From this expression we observe that the invariant goes always to zero as r — > oo: 

7 — >0(r" 6 ) , (20) 
as r — > oo 
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unless m = — 1 and n = 1. In this case, 1, 1) = and also 72 ( — 1 5 1) = 0. Thus, 

7(— 1,1) = 48r^r~ 8 . (20) 
The metric reads particularly simple in this case: 

ds 2 = -dt 2 + dr 2 + r 2 A(r)dQ 2 . 
We recall that for the Schwarzschild metric (m = n = 0): 

t 48r ° f9-n 

We are also interested in studying the behavior of the scalar invariants as r — > 2r . From 
expression (16) we see that 

I — 0[(r - 2ro)- 2 " 2 ™] , (22) 
as r — > 2r 

thus, for not occurring a singularity at r = 2r we must have: 

?7. < — 1, no singularity at r = 2r$ . (23) 

This condition can also be obtained from asking a non-singular behavior of the scalar 
curvature R, given by Eq. (14). Notice that n < — 1 makes singular at the horizon. 
However, this is only a coordinate singularity since the scalar invariants, as we have seen, 
are all finite on the horizon. 

One additional non-singular case is given when the term between curly brackets in 
Eq. (16) vanishes, i. e., 

m = n = 0; the Schwarzschild metric . (24) 

The other interesting value of the radial coordinate to study is r = 0. In this case, we see 
that: 

I^O(r- 6+2n ) , (25) 
as r — > 

except when 7i = 0, that is for n = —m = 2. In this case also I 2 = 0, then 

4-Rr 2 

7(-2,2)=-/ , (26) 



as r^0 
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which is the same as in the Schwarzschild case. In fact, it is easy to show that for n = 
—m = 2 the the metric (7) can be carried into Schwarzschild form. The transformation 
X = r — 2tq and the identification of — tq with the mass M make the job. 

We would like now to study the occurrence of an event horizon at r = 2tq. Let us 
first observe that the Killing vector £( t ) = ^/—good t = A(r)^^d t becomes null at r = 2r 
when m + 1 > 0. We can thus study the outgoing null geodesies from r > 2tq and see 
under which conditions r = 2r is an outgoing null surface. 

The first integral of the geodesies motion (related to the time and two angular vari- 
ables) in our spherically symmetric gravitational field can be written as: 

»= \ , (27) 

r 2 MrY% = J , (28) 

A(r) m+1 ^- =E , (29) 
clA 

^] =A(r)- n+1 {E 2 A(r)- m - 1 -J 2 r- 2 A(r)- n + e} , (30) 



where e = 0, ±1 for null, spacelike and timelike geodesies respectively. 

We can describe the radial part of motion in terms of the effective potential. Then, 
for null geodesies, we define the impact parameter b as: 

• (31) 

The critical impact parameters b c for which photons with b > b c can escape to infinity and 
with b < b c are absorbed by the black holes, can be found to be : 



dV f 



eff 



dr 



dr 

= ' Tx^ 



= 



d 2 V eff 



dr 2 



<0 . (32) 



and 



Thus, the radial coordinate of the critical periastrom r c is: 

v 

— = 3 + m — m , (33) 
ro 

bc = r f • (34) 

(1 + m — n) 2 

An observer at rest in our gravitational field measures the velocity of a photon relative to 
his orthonormal frame 

v + ~ V=g^dt/dX ~ r [r) ■ (35) 
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A photon at r < r c will eventually scape to infinity instead of being trapped by the black 
hole at r = 2r if v? is positive and: 



})q m — n-\-l 

sin 6 < — Air) 2 , 
r 

where 5 is the angle between the propagation direction and radial direction. 

Thus, we conclude that the surface r = 2ro will act as an event horizon whenever: 



m - n + 1 > 



(36) 



An alternative derivation of the horizon properties of the surface r = 2ro can be 
obtained by the study of the outgoing radial null geodesies. In fact, the time spent by a 
photon emitted at ri to reach 77 as measured by an observer at infinity is given by 

At= / dt= / A(r)^^dr = rA(r)^ +0(A(r)^ +1 ). 

in L J ri V / 

when m — n > as r\ — > 2ro the photon will need a At — > 00 to leave the horizon 
neighborhood, thus indicating the presence of an event horizon at r = 2r$. For n < — 1 
(eq(23)) we have not singularities on the surface r = 2tq. However, goO diverges there thus 
giving an infinite horizon area. This is only a purely geometrical divergence bringing no 
physical consequences. In fact, we have seen that the surface r = 2ro effectively acts as an 
event horizon with respect to null rays. For massive particles there neither any inconvenient 
to enter in to the black hole in a finite proper time since its effective potential, V e ff. 



E 



Mr) 



l-n 



E 2 A(r)-( m+1) - ^A(r)~ n - 1 



remains bounded at and outside the horizon. Besides, tidal effects on the horizon are finite 
since curvature tensors are well behaved there. 



3. Astrophysical applications and discussion 

For not only dealing with the mathematical aspects of the solution and to obtain 
further restrictions on m and n, we will briefly study some astrophysical consequences. In 
particular, we will compute some physical quantities and show how much different they 
are from the General Relativistic results. Thus, confirming that metric(2) is indeed not 
the Schwarzschild one. 

Nature has the final word to decide between mathematical models. Thus, to see 
if the family of solutions of Brans-Dicke equations could represent nature, we can start 
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by computing its post Newtonian parameters (PPN). For an static spherically symmetric 
metric we can write the PPN metric as I 12 ^ 



ds 2 = 



(M\ n (M 
2 — + 2/3 — 
\ r ) \ r 



dt 2 + 



l + 2 7 



M 
r 



(dr 2 + r 2 dn 2 ) , 



(37) 



where (3 and 7 are two of the ten PPN parameters measuring, respectively, the amount of 
nonlinearity in the superposition law for #00 and the amount of space curvature produced 
by the unit rest mass. 



By transforming our metric (7) to isotropic radial coordinates, r, 

r = f(l + f /f) 2 ; f = r /2 , 

we find W 



ds 2 = 



r /r 



— \ 2(m+l) 



1 + r /r 



dt 2 + (1 + f /f) z 



1 



- /-\ 2n 

r /r 



1 + r /r 



(rff 2 + f 2 d0 2 ) 



(38) 



(39) 



By expanding the coefficients of this metric and comparing them to those of Eq. (37) we 
obtain: 

1 — ra 

0=1 ; 7= -r ; M = (m + l)r . (40) 

m + 1 

Thus, when m — > —n we have agreement with the solar system experiments. In particular, 
results of time delay measurements gives t 10 ': I7 — 1| < 10 -3 . 

To find observational differences between metric (7) with m — > — n (see (1)), and the 
Schwarzschild one, we must, then, look at strong gravitational field effects. We study some 
of such effects in accretion disks, scattering of photons and Hawking radiation. 

The standard model of galactic hard X-ray sources is a binary stellar system formed 
by a normal star transferring matter onto its companion star, which is a compact object. 
This matter, falling inward in quasi-circular orbits, will form an accretion disk, which will 
emit the observed X-rays. 

The friction due to viscosity will generate heat, which is radiated away through the 
disk surfaces. This energy is supplied by the loss of the total energy of the gas, while 
going through the disk, down to the last stable circular orbit. After this, the gas would 
fall almost without radiating [ 13 L 

Using Schwarzschild' s metric (see refJ 14 ! also for the Kerr case), the last stable cir- 
cular orbit has a radial coordinate r c = 6M, where M is the black hole mass. At this r c 
the energy "at infinity" per rest energy is E c = (8/9) 1 / 2 . If we take E = 1 at the external 
radius of the disk and a steady flux of matter (or its temporal average) the total luminosity 
of the accretion disk will be: 

L=(1-E C )M , (41) 



with M = mass per unit time entering the disk. 
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When we compute E c from metric (7) we obtain (for — n large): 



tBD 

Jgr = °- 958 ( 42 ) 

Another potentially observable quantities are the orbital frequency of the last stable cir- 
cular orbit (as seen by an observer at infinity) and the redshift at infinity which are given 

by[15] 

BD Z BD 

^ = 0.931 ; -^ = 0.936 , n -> -oo (43) 

As we have seen, the metric (7) is in agreement with the solar system observations 
and experiments when m — > —n. In addition, n < — 1 for having a regular horizon. Thus, 
Eqs. (42)-(43) give results close to those of General Relativity. 

Unfortunately, present uncertainties in the modeling and observation of accretion 
disk do not provide accurate enough data to discriminate between metric (1) and the 
Schwarzschild one. 

From the study of null geodesies we made in the last section we can obtain the total 
scattering cross section for photons: 

, 2 M 2 (3 + m - n) m " n + 3 

a = nbi = n- — -j- { — . (44) 

(m + 1) 2 (1 + m - n) m - n+1 v ' 

Let us observe that m = n gives the same results as for a Schwarzschild black hole. This 
is so, because when m = n the metric (7) can be written conformal to the Schwarzschild 
one, i.e. ds 2 = A(r) n ds 2 Schw and light rays do not "feel" conformal factors. 

When we compare this cross section to the General Relativistic result, a GR = 277rM 2 , 
in the case m — > — n and n < —1, we find that, again, the results are close to those 
produced by the Schwarzschild metric. As n goes to more negative values we have a quick 
convergence to the asymptotic value: 



a GR 



1.095 , n -> -oo (45) 



From the results above one sees that the Brans-Dicke gravitational field studied 
seems to be weaker than the Schwarzschild one. This conclusion will be reinforced when 
we compute the surface gravity on the horizon. Here we find the relatively strongest 
difference from the Schwarzschild's results. The surface gravity plays an important role 
when one studies the thermodynamics of black holes because it is related to the temperature 
associated to quantum effects close to the horizon. For a static spherically symmetric 
system it is given by: 

K=-\ , 9 ' 00 =(m+l)^(r)^ . (46) 
2 v-googrr r z 
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When we evaluate it at r = 2r we find: 

{0, for m > n 

oo, for m < n , (47) 

! Sr = TM = k Schw, for m = n 

thus, we obtain the Schwarzschild value for the conformal case m = n. 

We can write also our metric in terms of Kruskal-like variables. To this end, let us 
define first the null variables u and v by 

du = dt-dr* ; dv = dt + dr* , (48) 

where: 

dr* =rA(r)^- 1 dr , (49) 

and then to 

u = — exp(— Khu) ; v = exp(KHv) , (50) 
where Kh is the surface gravity evaluated at r = 2r$. 
Finally, we obtain ^ 

ds 2 = -A(r) m+1 K~ 2 exp(-2K H r*)dudv + r 2 dQ 2 . (51) 

When r — > 2r$ and m = n we have 

^(2r )=e-+ 1 f^-'j . (52) 
\m+lj 

So, metric coefficient are finite on the horizon. 

When we bring together all the conditions for having a regular black hole, we obtain: 
n < — 1 for the horizon not being a singular surface (Eq. (23)) and m — n + l>0forr = 2r 
acting as an event horizon (Eq. (36)). If in addition, we ask that the solution should be 
in agreement with the observations carried out in the solar system, the PPN parameters 
should coincide with those of General Relativity with great precision. As we have seen, 
this is achieved when \u\ — > oo, i. e. m + n = 0. This, in turn, gives a constant scalar field 
outside the horizon (see Eq. (8)). It is notably that in this case, the Ricci tensor has one 
of its components different from zero, i. e. 

# n =2n(2-n)pi (l-2^)" 2 . (53) 

This is so, because in spite of the scalar field <p being constant and thus its derivatives going 
to zero, the coupling constant |u;| goes to infinity in such a way that the product appearing 
on the right hand side of the field equations (5) gives a finite value, i. e. (Eq.(53)). This 
fact has very important consequences for the Hawking theorem ^ establishing the identity 
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of Brans-Dicke and General Relativity Black Holes. Indeed, the limiting case m + n — > 0~ 
is contained within our family of solutions and it is well defined (For example m = —2, 
n = 2 gives the Schwarzschild metric). In this black hole solution (1): 

ds 2 = -A[rY- n db 2 + Air^dr 2 + r 2 A{r) n dQ 2 , 

the parameter n plays the role of a classical Brans-Dicke hair. It has its origin in the 
particular coupling of the Brans-Dicke scalar field. Their effects at large distances can 
be absorbed in a redefinition of the mass of the black hole and thus as we have seen, at 
Post-Newtonian level this metric coincides with the Schwarzschild one. However, as we 
study strong gravitational field effects, their results are dependent on the value of n. In 
some sense, thus, n has an intermediate range of action. 

When <p is not constant, black hole solutions are still possible due to the fact that the 
surface integral J ((p 2 ) ,a d'E a (where (p = <j) — <j)o) assumed to vanish t 3 ' 10 ' (under the implicit 
supposition of T^l^V > 0), here gives a non - zero contribution, i. e., 47rro</>o( m + n ) < 0' 
thus compensating the positive value of the integral J (</?,a) 2 \Z^gdx 4 . This is indeed so 
due to the particular form of the scalar field Eq. (8), which produces a stress tensor that 
violates the weak energy condition (with (m + n) < ensuring regularity of the field on 
the horizon). The finiteness of the surface integral can be understood by the fact that <7## 
for n < — 1 diverges on the horizon, thus producing a finite result when multiplied by the 
vanishing scalar field terms and integrated over the horizon surface. 

It is worth to stress that as m + n — > 0~ and n < —1, then u — > — oo. This 
is perfectly acceptable because there is no theoretical reason to restrict u> to positive 
values' 17 ' and experiments are consistent with |u;| > 500. For the allowed range of values 
of the parameters m and n (given by eqs. (23) and (36)), — oo < u < —4/3. Let us 
remember that the string theory selects the value' 2 ' us = —1; while the graceful exit 
problem is solved for W UJEI < 20. 

The no hair theorem can be overcome because the weak energy condition is violated 
by the energy momentum tensor of the Brans-Dicke field. In fact, 

Too = n(2-n)4(l-2r /r) m - n , (54) 

that for n < —1, takes always negative values (independent of the limit m — > — n). 

Another interesting result is that for the subfamily (1), in particular (in general see 
Eq. (36)), the surface gravity will be zero (see Eq. (47)). Hence, these black holes are 
truly "black", even at the semiclassical level, in the sense that not Hawking radiation is 
expected to take place here. It is worth to remark here that the divergence of the horizon 
surface not only does not affect the computation of relevant physical quantities, but can 
also be interpreted, with regards to the thermodynamics of black holes, as suggesting an 
infinite entropy for our black hole solutions. This, in turn, is consistent with its associated 
semiclassical zero temperature. 
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We remark that this kind of analysis can be also carried out for the generalization of 
the Kerr-Newmann metric in the Brans-Dicke theory I 11 !. 

The problem of stability of solution (1) is now under study by the present authors, 
but we can advance some comments: Matsuda ^ has found, studying the spherically 
gravitational collapse of a star in Brans-Dicke theory, that it does not necessarily produce 
a Schwarzschild black hole, but can also produce the black hole solution given by metric 
Eq. (1). 

The radiation of the scalar field will be damped by a factor t 3 ' 10 ! (2 + cj) _1 , which 
vanishes as |u;| goes to infinity. Thus, we think that metric (1) is a viable candidate to 
represents the black holes in nature. 
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